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fN| ' Measuring the primordial matter power spectrum is our primary means of probing 

' unknown physics in the very early universe. We allow the primordial power spectrum to 

J> ■ be an arbitrary function, and parametrize it in terms of wavelet band powers. Current 

cosmological data correspond to 11 such wavelet bands. We derive constraints on these 

band powers as well as Hq, f^/i 2 and il m h 2 from current Cosmic Microwave Background 

Anisotropy (CMB) data using the Markov Chain Monte Carlo (MCMC) technique. Our 

£f) [ results indicate a feature in the primordial power spectrum at 0.008 < ^/(/iMpc -1 ) < 

0.1. MAP and Planck data should allow us to put tighter constraints on the primordial 

' power spectrum. 

Q-f 

6 ■ 
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52 ' 1- Introduction 

. £h ! As a result of recent cosmological data, inflation (Guth 1981; Albrecht & Steinhardt 1982; 

Gott 1982; Linde 1983) has become increasingly well established as the plausible solution to the 
problems of standard cosmology (Kolb & Turner (1990); see Peebles & Ratra (2003) for a recent 
review). The primordial power spectrum is our primary window into unknown physics during 
inflation (Wang, Spergel, & Strauss 1999; Chung et al. 2000; Enqvist & Kurki-Suonio 2000; Lyth, 
Ungarelli, & Wands 2002). It is of critical importance that we try to extract the primordial power 
spectrum, Pi n (k), from observational data without assuming specific forms for it. 

Cosmological parameters are being measured to impressive precision with the help of recent 
CMB and large scale structure (LSS) data. The parameter constraints thus deduced are however 
sensitive to assumptions regarding the power spectrum of primordial density perturbations (Kinney 
(2001); Souradeep et al. (1998)). The primordial power spectrum is often assumed to be a power 
law, which represents many inflationary models (for example, see Linde (1983); Freese, Frieman, 
& Olinto (1990); La & Steinhardt (1991)). With such a parametrization the primordial power 
spectrum has been found to be scale invariant to a very good approximation, and its amplitude 
constrained (see for example Gorski et al. (1994); Lewis & Bridle (2002)). However, there are 
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many inflation models that predict primordial power spectra which cannot be parametrized by a 
simple power law (for example, Holman et al. (1991ab); Wang (1994); Randall, Soljacic, & Guth 
(1996); Adams, Ross, &; Sarkar (1997); Lesgourgues, Polarski, & Starobinsky (1997); Lesgourgues 
(2000)). These can represent unusual physics in the very early universe. For example, inflation 
might occur in multiple stages in effective theories with two scalar fields (Holman et al. 1991ab), or 
in a succession of short bursts due to symmetry breaking during an era of inflation in supergravity 
models (Adams, Ross, & Sarkar 1997). 

With the quality of data improving, more attention is being paid to the nature of the primordial 
perturbations (for example, see Covi, Lyth, & Melchiorri (2002); Leach & Liddle (2002)). As 
more observational data become available, they increase our ability to probe the primordial power 
spectrum, Pj n (/c), as an arbitrary function of scale. A model independent determination of the 
Pin(k) could uniquely probe physics of the very early universe, test what we have assumed about 
early universe physics, and provide powerful constraints on inflationary models. Wang, Spergel, & 
Strauss (1999) explored how this can be done with the CMB data from MAP 1 and LSS data from 
SDSS, using a piecewise constant function for Pi n (k). Wang & Mathews (2002) used the CMB data 
from Boomerang, Maxima, and DASI to place constraints on Pi n {k) using linear interpolation to 
approximate the function between several k values equally spaced in logfc. 

Here we employ wavelets in a model independent parametrization of the primordial power 
spectrum and obtain constraints from current CMB data. We use only CMB data here, but LSS 
(for example, see Hamilton & Tegmark (2002); Percival et al. (2002); Bahcall et al. (2003)), and 
CMB polarization (for example, see Kovac et al. (2002)) data can all be added to help break or 
reduce degeneracies between different cosmological parameters and help better constrain Pi n {k). 

We describe the wavelet parametrization of the primordial power spectrum in Sec. 2 (see Ap- 
pendix A for further discussion). In Sec. 3, we discuss the techniques that we have used to optimize 
our method. Sec. 4 contains our results. Sec. 5 contains a summary and discussions. 

2. Wavelet Band Powers of the Primordial Power Spectrum 

Assuming that the primordial matter density fluctuations form a homogeneous Gaussian ran- 
dom field, the wavelet band powers of the primordial power spectrum Pj„(/c), can be written as 
(see Eq.(AlO)) 

2 

Pinik), (1) 

where ip is the Fourier transform of the wavelet tp (see Appendix A). Pj's represent a scale-by-scale 
band averaged Fourier power spectrum, where the banding is not arbitrary but well defined, and 
naturally adaptive (as wavelets by construction afford better k resolution at smaller k). The Pj's 
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are mutually uncorrelated by construction (see Appendix A and Fang & Feng (2000) for further 
discussion of these issues). 

To illustrate, Figure 1 shows a primordial power spectrum with features (solid line), similar to 
that discussed by Elgaroy, Gramann, & Lahav (2002). The points are the wavelet band powers Pj, 
computed using Eq.(l). These Pj's can be thought of as being measured exactly from ideal data. 
Clearly, the wavelet band powers of the primordial power spectrum are excellent approximations of 

2 

the primordial power spectrum at the central k values of the wavelet window functions, i/j (^-) 
(dotted curves in Fig.l). Since the window functions are the modulus squared of the Fourier 
transform of dilations, by factors of 2, of the basic wavelet (see Appendix A), the peaks of these 
window functions are separated by factors of 2 in k. While we can add more bands to lower k, and 
3 more bands at higher k that go up to 2tt in k, and reconstruct the amplitude of the power in 
these additional bands, the locations of the bands or their widths cannot be changed. This is why 
in this scheme the banding is not arbitrary. 



Hence we parametrize the primordial power spectrum Pi n {k) as follows [see Eq.(All)] 



(2) 



The wavelet band powers Pj can be estimated from data that depend on the primordial power 
spectrum P in (k), for example, CMB and LSS data. Here we use only CMB data and estimate 
cosmological parameters together with the wavelet band powers Pj , which represent a measurement 
of the primordial power spectrum at the central k values of the wavelet window functions. P,'s 
of unity corresponds to the default power of A 2 S = 2 x 10 -9 , and if all P^-'s are set to unity we 
recover exactly the C/ spectrum that results from a scale invariant primordial power spectrum 

*.(*) = 4 

2 

Figure 2 shows how each wavelet band power window function, ip (i) [shown in Figure 1], 
maps on to a window function in the CMB angular power spectrum multipole number I space, for 
one set of cosmological parameters. These window functions have been numbered 7 to 17, which 
correspond to a range of 0.0001 < /c/(/iMpc -1 ) < 0.5. The 11 band powers, P-j through Pyj (from 
small to large k), are sufficient to constrain the primordial power spectrum from current and near 
future cosmological data. CMB data, up to an l max of 1500, are insensitive to band powers Pi to 
Pq, which correspond to smaller k and have been set equal to P?, and to band powers higher than 
Pi7, which correspond to larger k and have been set equal to Pyj. The solid curve is the CMB 
angular power spectrum C\ that includes contributions from all the Pj's [set to unity in this Figure]. 
Given that the CMB data are in the form of band powers in C[S, and each multipole I maps on 
to a range in k (for example, see Figure 4 of Tegmark, & Zaldarriaga (2002)) with dependence on 
cosmological parameters, we see that in constraining Pj using cosmological data some correlations 
between the P,'s will be inevitably introduced. The P,'s are also correlated somewhat with the 
cosmological parameters in ways that can be understood from Figure 2. 
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For a comparison of different banding methods, instead of wavelet band powers we can parametrize 
the primordial power spectrum as a continuous and arbitrary function determined by its amplitude 
at several wavenumbers, equally spaced in log(fc), using linear interpolation to approximate the 
function at intermediate wavenumbers. This is the linear interpolation binning method of Wang &i 
Mathews (2002). In this work, we determine amplitudes at 11 wavenumbers corresponding to the 
central k values of the wavelet band powers. Figure 3 shows how each of these "binned" amplitudes 
map onto the CMB angular power spectrum multipole number I. The solid curve is the CMB 
angular power spectrum C\ with contributions from all the bins [all bin amplitudes set to 1]. The 
power in each bin, denoted here by aj, is correlated with the power in neighboring bins from the 
start, and ultimately the resulting correlations will be different from the wavelet banding method. 
The error bars on the estimates will also be correlated. See Wang & Mathews (2002) for further 
details about this method. 



3.1. Optimization Techniques: Wavelet Projections and Markov Chain Monte Carlo 

The usual approach to deriving cosmological constraints from CMB data is to grid over all the 
parameters being estimated and compute the likelihood at each point. Here, besides the 11 wavelet 
band powers (Pj's), we constrain Ho, Qbh 2 and Q c h 2 , assuming a fiat universe with a cosmological 
constant A, and ignoring reionization and tensor modes. This results in 14 parameters. In order to 
place the most accurate and reliable constraint on the primordial power spectrum, we have chosen 
to compute the theoretical CMB angular power spectra at the accuracy of CMBFAST (Seljak & 
Zaldarriaga 1996) and CAMB (Lewis, Challinor, & Lasenby 2000). Our analysis would have been 
extremely time-consuming even on supercomputers. Fortunately, we are able to use two techniques 
which have made it possible for the analysis to be completed in a timely fashion. 

First is the wavelet projections of the CMB angular power spectrum, C/'s. Using Eq. (2), we 
expand the CMB angular power spectrum as follows: 



where the cosmological model dependent transfer function A<n(/c,r = To) is an integral over con- 
formal time r of the sources which generate CMB fluctuations, To being the conformal time today, 
and s represents the cosmological parameters other than the Pj's. We use CAMB 2 to compute 
the CMB angular power spectrum, in a form such that for given cosmological parameters other 
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Fig. 1. — An example of a primordial power spectrum (solid line), and its wavelet band powers Pj's 
(-7=7,17; points), together with the corresponding wavelet window functions (dotted lines). The y 
axis is in arbitrary units. 
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Fig. 2. — Mapping of the wavelet window functions of Figure 1 into window functions in the CMB 
multipole I space (dotted curves). The solid curve is the C\ spectrum that is the sum of contributions 
from all the wavelet bands (all the band powers, Pj, are set to unity here for illustration). 
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than the Pj's, the C/(s) are computed, so that there is no need to call CAMB when we vary the 
Pj's. This results in significant computational speed up. We expect this technique to be essential 
in constraining the primordial power spectrum with sufficient detail and speed. 

The second technique is the Markov Chain Monte Carlo (MCMC) technique. The large number 
of parameters being varied here necessitates the use of this technique in the likelihood analysis. At 
its best, the MCMC method scales approximately linearly with the number of parameters. The 
method samples from the full posterior distribution of the parameters, and from these samples 
the marginalized posterior distributions of the parameters can be estimated. See Neil (1993) for 
a review, and Hannestad (2000); Knox et al. (2001); Lewis & Bridle (2002); Rubino-Martin et al. 
(2002) for applications of this method to CMB analysis. 



3.2. Recipe for the Analysis 

The steps followed in the analysis are as follows: 

(i) Establish the multi-dimensional parameter space that span all the parameters to be estimated 
from data. In our case, the total number of parameters is 14, including the wavelet band powers 
P/s. 

(ii) Start at some point in our 14d parameter space, which corresponds to an initial set of param- 
eters. Note that the MCMC chain soon loses memory of this starting point and it can be verified 
that the result is independent of it. 

(iii) Find the likelihood of this set of parameter values given the data: 

-2ln{L) =x 2 = Yl [ /n (Cfc + Xb) ~ HC l b h + x b j\ G h , v [ln(C v + x v ) - ln(C$ + x v j\ , (4) 
b,b' 

where C b and Cf 1 are the experimental and theoretically estimated band powers respectively, G^y 
is the corresponding band-band correlation matrix and x b are the offset parameters in the offset- 
lognormal ansatz of Bond et al. (1998), or for experimental results computed without this 
ansatz. For this purpose, we calculate the C\ spectrum for the given set of parameters, estimate 
the expected experimental band powers (for the given set of parameters) using the band power 
window functions available for the different CMB experiments, and compare these estimates to 
the measured band powers, taking into account the instrumental noise. We use offset-lognormal 
band-powers whenever available, and analytically marginalize over known beam and calibration 
uncertainties for each experiment, as described in Bridle et al. (2002). 

(iv) Take a random step in parameter space. The probability distribution of the step is taken to 
be Gaussian in each direction with an rms that is neither too small (or the chain will exhibit poor 
mixing) nor too large (or the chain efficiency will be poor). Ensure that the new point lies within 
the wide priors that define the likelihood region that one chooses to explore (our priors are defined 
in the next section). Calculate the likelihood at this new point (which corresponds to a new set of 
parameters). 
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(v) Use the Metropolis-Hastings algorithm (an MCMC sampler based on the Metropolis-Hastings 
algorithm has been made available in the software package CosmoMC (Lewis & Bridle 2002)) to 
determine whether to accept this point or to discard it and propose another point in the parameter 
space. 

(vi) Repeat steps (hi) to (v) until a large number of samples from the posterior distribution of these 
parameters have been chosen. 

(vii) Find the Id marginalized parameter distributions and confidence limits for each parameter. 

(viii) Check for convergence by running several chains starting at different initial values giving the 
same final distributions, and a couple of longer chains to ensure good sampling. 

It would be interesting to include reionization and curvature in the analysis, and obtain ad- 
ditional constraints from different kinds of LSS data and ultimately also CMB polarization. We 
defer such extensions to the near future. 

4. Results and Discussion 

We use most of the recent high precision data for which experimental details are publicly 
available. These are the latest Boomerang (Ruhl et al. 2002), Maxima (Lee et al. 2001), DASI 
(Halverson et al. 2002), VSA (Scott et al. 2002), CBI (Pearson et al. 2002), ACBAR (Kuo et al. 
2002) and Archeops (Benoit et al. 2002), together with COBE (Smoot et al. 1992). We estimate 
the wavelet band powers of the primordial power spectrum, as well as Hq, Q^h 2 and Cl c h 2 . 

All results shown in this paper are for the wavelet Daubachies 20 [see Appendix A]. 

Allowing all 14 parameters to vary, the results obtained are shown in Figure 4. We have 
used a weak prior of the age of the universe to > 10 Gyrs. All the Pj's are consistent with unity 
at ~ 1.5o". However, there is some indication of a feature in P{ n (^k^ though at low significance. 
Figure 4 shows that current CMB data seem to favor a dip at a k of ~ 0.01 /iMpc -1 (from P12) 
and excess power k of ~ 0.02 to 0.03 /iMpc -1 (from P13 and P14), followed by another small dip 
from Pi6 at k of ~ 0.1 /iMpc -1 . The cosmological parameters are constrained to h = 0.56 ± 0.09, 
n b h 2 = 0.020±0.005 and Q c h 2 = 0.161±0.028. 3 , 4 Note that the derived value of the Hubble constant 



3 In this paper we quote the mean of the derived Id distributions, instead of the maximum, for the constrained 
parameters, following Lewis & Bridle (2002). As discussed in Lewis & Bridle (2002), the MCMC samples from the 
posterior do not provide accurate estimates of parameter best-fit values, because in higher dimensions the best-fit 
region typically has a much higher likelihood than the mean but occupies a minuscule fraction of parameter space. 
Our main results are not affected by this. 

4 This corresponds to fl b = 0.062(0.048,0.079] and tt c = 0.51[0.42, 0.60]. Note the degeneracies between the 
cosmological parameters considered here. Due to the geometrical degeneracy, the location of the first acoustic peak 
nails the curvature of the universe (here taken to be zero) but a degeneracy between fi m and $7a, or equivalently O m 
and h remains. The height of the first peak is degenerate for Q,bh 2 and fi c /i 2 . Precise determination of the second 
and third peak heights can help determine Q b h 2 and thus ease degeneracies (for example, see Efstathiou & Bond 
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h is lower than the value of h derived assuming a scale-invariant primordial power spectrum (i.e., 
setting all the Pj's to the same constant value). Assuming scale-invariance gives h = 0.71 ± 0.06, 
n b h 2 = 0.023±0.001, Q c h 2 = 0.122±0.016 and all P j , s= 0.90±0.06. A chi square analysis indicates 
that a model with the P,'s estimated from data (i.e., allowing Pj n (fc) to be a free function) and the 
best fit scale invariant model (setting all P,'s equal to a constant) fare about the same. For a low 
h, the derived Q m is high. If we used an external constraint, such as from LSS, to keep Vt m low, 
the derived h would be higher; this is a result of parameter degeneracies in CMB data. 

Figure 5 shows the CMB power spectrum corresponding to the fitted cosmological parameters 
and Pj's and also the CMB power spectrum for the same cosmological parameters but with all P,'s 
set to unity. One can see where the Pj's are contributing. Most distinctly, by allowing the P,'s to 
vary, the CMB angular power spectrum receives significant positive contribution around the first 
peak. 

The corresponding results for the case of the linear interpolation binning method are in Figure 
6. Note that the primordial power spectrum reconstructed in this way is consistent with the 
primordial power spectrum reconstructed using wavelet band powers [Figure 4]. This is reassuring 
since in the two methods the parameters being reconstructed have different correlations amongst 
themselves and with the other cosmological parameters [see Figures 2 and 3]. The bin amplitudes 
a\ to clq do not rule out the excess power at small k that was found by Wang & Mathews (2002). 
But from the marginalized distributions we see that these bin amplitudes are not constrained and 
vary over the entire range that was allowed in the runs. In comparison the wavelet banding method 
works better. This is as expected since the wavelet band powers are uncorrelated by construction, 
hence are less correlated when estimated from data than the linear interpolation bins. 

Now we examine the effect, on the wavelet band power results, of imposing statistical priors. 
If a Gaussian prior is introduced in the value of h, a weak conservative prior of 0.6 ±0.1 (Branch 
1998), together with the above mentioned weak age prior, what results is almost identical to Figure 
4. If a stronger prior on h is introduced, that of h = 0.72 ± 0.08 (Freedman et al. 2001) we find 
that the feature in the 0.008 < A;/(/iMpc _1 ) < 0.1 range remains, though it reduces in amplitude 
as a whole by nearly 10% (Figure 7). Cosmological parameters are simultaneously constrained to 
be h = 0.64 ± 0.06, n b h 2 = 0.023 ± 0.004 and Q c h 2 = 0.147 ± 0.022. Note that when the P/s were 
varied near the currently favored values of the cosmological parameters (h = 0.68, Q b h 2 = 0.024, 
and Q c h 2 = 0.12), the same form for the primordial power spectrum, with the whole spectrum 
shifted down by about 20%, was obtained. Without priors on the Hubble constant, the current 
CMB data seem to prefer a somewhat low Hubble constant of h = 0.56 ± 0.09 with the freedom 
that is allowed in the shape of the primordial power spectrum. It is not clear whether this is due to 
any systematic effects. Note that however, this lowish value of the Hubble constant is within the 
uncertainties of other independent measurements of Hq (for example, see Branch (1998); Freedman 
et al. (2001); Gott et al. (2001); Saha et al. (2001)). 



(1999); Hu et al. (2001); White & Colin (2002)). 
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We note that several other groups have studied the constraint on specific types of features of 
the primordial power spectrum from current observational data; our findings are consistent with 
these results (Silk, & Gawiser 2000; Atrio-Barandela et al. 2001; Barriga et al. 2001; Hannestad, 
Hansen, &; Villante 2001; Elgaroy, Gramann, & Lahav 2002). 

The advantage of our approach in this paper versus previous work is as follows. We parametrize 
the primordial power spectrum by its wavelet band powers Pj. This is a model independent 
parametrization of Pi n (k) 5 . Wang, Spergel, & Strauss (1999) parametrize the primordial power 
spectrum by a step function [i.e., top-hat banding]; this leads to an estimated primordial power 
spectrum Pi n {k) that is discontinuous, while the Pi n {k) estimated using wavelet band powers is 
continuous. Wang & Mathews (2002) used linear interpolation of Pi n {k) values at several k values 
[equally spaced in logfc]; this lead to a continuous estimated Pi n (k), but the Pi n {k) values estimated 
from data are strongly correlated. The wavelet band powers Pj are mutually uncorrelated by 
construction [see Eqs.(l)(2), though some correlations are inevitably introduced because the data 
that are used to constrain them are bands in I space], and the banding in this method is not 
arbitrary. Note also that by computing the wavelet projections of C\ [see Eq.(3)] for each set of 
cosmological parameters excluding the Pj's, we avoid computing Ci when we vary the wavelet band 
powers Pj in the likelihood analysis. This, together with the MCMC technique, made it possible 
for us to estimate all the 11 relevant wavelet band powers from the current CMB data in a timely 
fashion. We expect that our method will be efficient in yielding tighter and detailed constraints on 
Pin(k) when applied to MAP and Planck data. 

Finally we note that the primordial power spectrum results presented here are for k in units 
of h Mpc -1 . These units are found to yield more constraining results and this is important for 
pre-WMAP data. Some published results are for Pi n (k) with k in units of Mpc -1 . Note that when 
Pin{k) is given with k in units of Mpc -1 , the feature we have found in this paper would be shifted 
to smaller k, to near k ~ 0.01 Mpc -1 . 

5. Conclusions 

A model independent determination of the Pi n (k) could uniquely constrain unknown physics 
in the very early universe, test what we have assumed about early universe physics, and provide 
powerful constraints on inflationary models. Thus rather than assuming specific forms for Pi n (k), 
we have used a wavelet band power expansion to extract Pi n (k) as a free function, using recent high 
precision CMB data. The wavelet band powers parametrization of the primordial power spectrum 
has the following features: in this scheme the banding is not arbitrary but well defined and adaptive. 
In terms of these band powers the primordial power spectrum can be reconstructed as a smooth 
function. The band powers are mutually uncorrelated by construction, and they are excellent 



5 Miller et al. (2002) do a model independent analysis of the CMB angular power spectrum C;'s. 
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approximations of the primordial power spectrum at the central k value of the wavelet window 
functions. Although in estimating these from CMB data, which are band powers in multipole 
space, some correlations are inevitably introduced because of the cosmological model dependent 
nonlinear mapping between wavenumber k and multipole I spaces. 

The wavelet band powers of Pi n {k) that we have extracted from current CMB data seem to 
indicate a feature in the primordial power spectrum at 0.008 < /c/(/iMpc -1 ) < 0.1, though only 
at low significance. A chi square analysis indicates that a model with the estimated Pj's and 
cosmological parameters and the best fit scale invariant model fare about the same. Future data 
will better help distinguish between these models. The linear interpolation binning approach of 
Wang h Mathews (2002) yields an estimated Pi n {k) with a similar feature at roughly the same 
location in k with comparable significance. Our results are consistent with previous work by Wang 
Sz Mathews (2002). MAP and Planck 6 data should allow us to put a tighter constraint on the 
primordial power spectrum (Lesgourgues, Prunet, & Polarski 1999; Wang, Spergel, & Strauss 1999; 
Hannestad 2001; Matsumiya, Misao Sasaki, & Yokoyama 2002; Tegmark, & Zaldarriaga 2002). 

We acknowledge the use of CAMB and CosmoMC. This work was supported by NSF CAREER 
grant AST-0094335 at the Univ. of Oklahoma. We thank the referee for helpful comments. 
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A. Wavelet Band Powers 

Here we describe our wavelet parametrization of the primordial power spectrum in detail. 

The wavelet transform bases are obtained from dilations and translations of a certain (mother) 
function ip(x) via 

M*) = (xj ^( 2 ^/ L -0- (Al) 

where ip(x) is in general real, defined on the interval [0,1], and obeys several restrictive mathe- 
matical relations first derived by Daubechies (1992) in order for the resulting wavelet basis to be 
discrete, orthogonal and compactly-supported. These are the kind of wavelets we consider here. 
See, for example, Press et al. (1994) for an introduction to wavelets, and Barreiro et al. (2000) 
and Tenorio et al. (1999) for applications of spherical wavelets to CMB data on the sky. The j and 
I are scale and position indices respectively, and the wavelet bases are orthogonal with respect to 
both these indices, 



oc 



tp ji l(x)tp j/i i / (x)dx = Sjj'Su'. 



A periodic function f(x) of period L, sampled at N = 2 J equally spaced points between and L, 
can be expanded in terms of the wavelet basis as 

J-l 2^-1 

f(xi) = E E (A2) 

3=0 1=0 

where the coefficients bjj are given by 

b jt i= [ f(x)^ tl (x)dx. (A3) 
J o 

The scale index j increases from to J — 1, and the wavelets with increasing j represent the 
structure in the function on increasingly smaller scales, with each scale a factor of 2 finer than the 
previous one. The index I (which runs from to 2 3 — 1) denotes the position of the wavelet ipj i 
within the jth scale. Thus bjj measures the signal in f{x) on scale L/2- 7 , and centered at position 
IL/23 in physical space and centered at wavenumber 2tt x 2 j, /L in Fourier space. 

The Fourier decomposition of function /(xj) is given by 

N-l 



f( Xi ) = E e n e i2nnx > /L (A4) 



n=0 



and the Fourier coefficient e n is 



e n = j£ f(x)e~^ L dx. (A5) 
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Since both the discrete wavelet transform (DWT) and the Fourier transform (FT) bases are com- 
plete, there exists a relationship between the Fourier and wavelet coefficients. Substituting equation 
(A4) in (A3) gives 

oo 

hi = S e «fe( _n )' ( A6 ) 

n=— oo 

and similarly equations (A5) and (A2) give 

j oo 2^-1 
3=0 1=0 

where ^j,i{n) is the FT of the wavelet ipjj and is related to the FT of the basic wavelet (using (Al)) 

by ' _ i2 

feW = £ Mx)e-*™*' L dx = (|) $ g) e-^l*. (A8) 
From (A6), the covariance in wavelet space is given by 

oo 

(bj,lbj>,l>}= ( e n> e n')fe( n )^]',i'( n ')- ( A9 ) 



n.n =— oo 



For a homogeneous Gaussian random field, ignoring the often very small correlations that may 
exist between the wavelet coefficients (see Frazier, Jawerth & Weiss (1991), Walter (1992), Zhang h 
Walter (1994), and Tenorio, Stark & Lineweaver (1999); such correlations, ignored in mosts works, 
reduce further with the regularity of the wavelet), (bjibj' i') = PjSj,j'8i,i' - 7 Also for a Gaussian 
random field the Fourier amplitudes (|e n |) have a Gaussian one point distribution and their phases 
are random, so that (e n ,e n /) = P(n)S n>n i. Thus 

1 °° 2 

p i = 2i £ p{n) - (A10) 

n=— oo 

where Pj is the variance of b^j, the power of perturbations in wavelet coefficients of scale j. 

The Pj's are thus the scale-by-scale band-averaged Fourier power spectrum, from which one 
can attempt to reconstruct the Fourier power spectrum as a smooth function. The weak correlation 
of the wavelet coefficients discussed above leads to an excellent approximation to P{n): 

oo 2 

P(n) = J2 P j\^{^)\ ■ ( A11 ) 

3=0 



7 When the random field is ergodic, the 2 J coefficients at a given scale can be taken as 2 J independent measurements. 
The average over 1 is thus a fair estimation of the ensemble average. 



Moreover, for a Gaussian random field the Pj's are (very nearly, see discussion above) uncorrelated: 

(^^■ + l) _ 2J+1 p^/2^ 

Scale-scale correlations, as defined above for order 2 (the value of the exponent), has been discussed 
in detail by Pando, Valls-Gabaud & Fang (1998), and Mukherjee, Hobson & Lasenby (2000), and 
others (note that there are half the number of coefficients at scale j than at scale j + 1). The 
existence of detectable scale scale correlations is an indication of mode-mode coupling and hence 
non-Gaussianity. 

In order to parametrize the primordial power spectrum in a model independent way, we need 
to estimate the power in certain bands in k. The bands should be logarithmically spaced in k, 
as in Wang & Mathews (2002). To find a unique number for the appropriate number of bands to 
use, we adopt the DWT approach to banding. While in the Fourier approach the phase space is 
split such that the resolution in wavenumber k is highest at all k (Ak — » 0), and the resolution in 
position x is lowest, (Ax — > oo), in the wavelet approach these resolutions are adaptive. We choose 
Ax oc 1/k, and Ak/k = log 2, so that an optimal chopping of the phase space is achieved whilst 
satisfying the uncertainty relation AxAk > 2tt (see Fang & Feng (2000) for a discussion). Wavelets 
afford good k resolution at small k and poorer resolution by factors of 2 as j (or position resolution) 
increases. Further, the wavelet band powers Pj, for a Gaussian random field, are uncorrelated by 
definition, and one cannot have more independent bands (Fang & Feng 2000). Equations (A10) 
and (All) show how the primordial Fourier power spectrum can be parametrized in terms of Pj's 
which represent a scale- by-scale band-averaged Fourier power spectrum with logio(2) spacing. 

Note the Eq.(AlO) is essentially just a suitable form of banding with window functions shown 
in Fig. 1. The equations preceding it in this section show that if the field in question, here a 
statistically homogeneous and isotropic Gaussian random primordial density field, were available 
to us then the Pj's would be the variance of wavelet coefficients of scale j. Although we have given 
equations for a Id field, if the primordial density field is isotropic the equivalence can be made. 

All results shown in this paper are for the wavelet Daubachies 20. We have also studied the 
case for the Symmlet 8 wavelet, obtaining very similar results. The larger the number associated 
with the wavelet, the more smooth is the wavelet in real space, and the lesser their compact support 
in real space (less localized, though compact support is technically more involved a concept than 
localization, and while wavelets can be localized in both real and Fourier space, it is impossible for 
a function to have compact support in both spaces). Since we hope to be able to pick up sharp 
features in the primordial (Fourier) power spectrum, wavelets that are smooth in real space are 
preferable for our purpose. In fact wavelets that have compact support in Fourier space rather than 
in real space (often called band-limited wavelets) should do better. Examples of such wavelets are 
the Shannon wavelet and the Meyer wavelet. We defer their use for a future paper as these are less 
frequently used wavelets so that the relevant software is not easily available. 
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Fig. 3. — Mapping of an arbitrary number of bins, here chosen to be 11, to correspond exactly to 
the central k values of the 11 wavelet bands discussed above, into window functions in the CMB 
multipole I space (dotted curves). The solid line is the C\ spectrum that includes contributions 
from all the bins (all the bin amplitudes are set to unity here for illustration). 
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Fig. 4. — The Id marginalized posterior distributions obtained upon varying all the 14 parameters. 
The bottom plot shows the constrained P^'s (j=7,17) versus scale. The dotted line indicates the 
best-fit scale-invariant model. The cosmological parameters are simultaneously constrained to be 
h = 0.56 ± 0.09, Q b h 2 = 0.020 ± 0.005 and Q c h 2 = 0.161 ± 0.028. 
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Fig. 5. — The solid line shows the CMB power spectrum at fitted values of the cosmological 
parameters and Pj's. The dashed line shows the CMB power spectrum for the same cosmological 
parameters at Pj's of unity. All the data points considered in the analysis are plotted. The error 
bars do not include calibration and beamwidth uncertainties. 
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Fig. 6. — Results obtained by parameterizing the primordial power spectrum with bin amplitudes 
ajs in 11 bins ( j = 1 , 1 1 ) corresponding exactly to the central k values of the wavelet bands. The 
dotted line indicates the best-fit scale-invariant model. The amplitudes a\ to a§ are unconstrained; 
hence their standard deviation is not meaningful. The cosmological parameters are simultaneously 
constrained to be h = 0.63 ± 0.10, n b h 2 = 0.021 ± 0.005 and Q c h 2 = 0.144 ± 0.029. 
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Fig. 7. — Same as Figure 4, but with a strong prior on Hq (see text). The cosmological parameters 
are simultaneously constrained to be h = 0.64±0.62, Q b h 2 = 0.021±0.004 and Q c h 2 = 0.144±0.022. 



